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We reduce the question whether a given quantum mixed state is separable or entangled to the 
problem of existence of a certain full family of commuting normal matrices whose matrix elements 
are partially determined by components of the pure states constituting a decomposition of the 
£ — i considered mixture. The method reproduces many known entanglement and/or separability criteria, 

and provides yet another geometrical characterization of mixed separable states. 
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!>■ ! I. INTRODUCTION 

Entanglement and separability problem. Entanglement is the most important quantum phenomenon, respon- 
i-£h sible for genuine, distinct and unique properties of the quantum world, and possibilities this world offers for future 
technological applications, such as quantum engineering, and quantum information Despite enormous efforts, 
many fundamental questions concerning entanglement remain open (for an excellent recent review see 1, and 1 
for some general geometric settings of the problem). In the seminal paper in 1989 Werner [J gave the definition of 
separable (i.e. non-entangled) states: a state of a bi-partite system is separable iff it is a mixture of pure product 
Q" 1 states. A simple question: given a state, is it separable or not, is known as the separability problem. Only in very rare 
l—— '■ instances we know operational sufficient and necessary criteria (SNC) that allow to solve this problem: 

• for 2 x 2 (two qubit) and 2x3 (qubit-qutrit) systems the SNC are given by the positive definiteness of the partial 



transform [5fl; this is the famous PPT criterion, introduced by Peres as necessary for separability in Ref. 



' • for 3 qubit symmetric ("bosonic") states PPT criterion is also SNC Q- 

• for continuous variables lxl (one mode per party) Gaussian states, PPT criterion (formulated at the level of 
correlation matrices) is a SNC [9[ . 

• for continuous variables mx n (all bipartite) Gaussian states there exist an operational SNC based on recursion 
for correlations matrices [To| . 

, ^ ' • for continuous variables tripartite lxlxl Gaussian states there exist an operational SNC based on " iteration" 

. of PPT condition for correlations matrices [111 ]. 

$H ' 

In general we have to rely either on only necessary criteria, or only sufficient ones, or on numerical approaches. 
Although there exist very efficient numerical procedures that employ optimization methods of semi-definite program- 
ming [12j, the complexity of the problem grows with the dimensionality of the underlying Hilbcrt spaces: in fact it 
has been proven that the problem belongs to the complexity NP-class [13| . 

Reformulations of the separability problem. The market for only necessary, or only sufficient criteria is 
growing constantly, and it is impossible to review it in a non-review style article (for this reason we recommend 
the readers the review 0)- There are also many attempts to reformulate the problem of separability in different 
mathematical terms. A paradigm example for such an approach is the formulation of the separability problem in 
terms of positive maps due to Horodeccy (BJ. A state is entangled iff there exist a positive map acting on, say, 
Alice, such that when applied to the state in question, it produces a non-positive definite operator. Similar approach 
deals with entanglement witnesses, i.e. observables that have positive averages on all separable states, but a negative 
average on some entangled state: A state is entangled iff there exist a witness operator that detects it, i.e. has a 
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negative average. Obviously both these approaches are not operational, but nevertheless they are extremely useful, 
since they allow to generate many necessary separability (sufficient entanglement) criteria via explicit construction of 
positive maps [l4| . witnesses [IB], and methods of their (local) measurements (cf. [lf|, respectively). 

We have recently presented another example of the reformulation of the separability problem employing harmonic 
analysis on compact groups [171 ]. In this approach quantum mechanical states are replaced by non-commutative 
characteristic functions defined on the considered group, and separability problem reduces to the question whether 
a characteristic function defined on a product group of two groups can be represented as a mixture of products of 
characteristic functions on each of the individual groups. The present paper is in a sense similar to the Ref. [l7l ]: 
we present yet another reformulation of the separability problem and reduce it to an apparently unrelated question 
of existence of a full (in a sense specified below) family of commuting normal matrices, whose matrix elements are 
partially determined by components of the pure states that constitute a decomposition of the considered mixed state. 

Decompositions of mixed states. A given genuine (not pure) mixed state p has infinite number of decomposi- 
tions in terms of projectors onto pure states. This fact has been already recognized by Schrodinger in 1935 [181 ] . and 
elaborated thoroughly from the more modern view by Hughston et al. [191 ] . Any decomposition of a density matrix of 
rank r into K projectors can be described in terms of a rectangular K x r matrix, whose r columns are orthonormal. 
Such objects are known in geometry to form a so called Vk,t = U(K)/U(K — r) Stiefel manifold [20| . Separability 
problem might be also formulated as a problem of statistical mechanics of a fictitious systems on the Stiefel manifold, 
characterized by a positive definite Hamiltonian (cost- function) that vanishes for separable states [21| . Here, we follow 
another avenue: we consider all decomposition of p into K terms for sufficiently large K: such decompositions are 
related via unitary transformations U{K). The matrix elements of the density matrix, on the other hand, form a 
Gram matrix of scalar products of certain vectors from this K dimensional space. The first chapter of the paper is 
thus devoted to the study of such Gram decompositions. It provides complementary results to the Ref. [l9j . 

Plan of the paper. As stated above, Section II is devoted to the Gram decompositions, and its main result is 
the Theorem 1, that describes how the two different Gram decomposition are connected. We present a reformulation 
of the separability problem in Section III, in the Theorem 2. Here, an example of so called Werner matrices [H is 
elaborated in detail. The section IV contains the main result of this paper: a novel (but unfortunately not immediately 
operational) SNC for separability in terms of existence of what we call a full family of commuting normal matrices 
(FFCNM) (Theorem 3). We specify this result to the particularly simple case of 2 x N systems, where the separability 
SNC requires existence of a single normal matrix, whose matrix elements are partially known. Here we use the general 
properties of the density matrices in 2 x N systems (as presented in Appendix) and formulate elegant theorems on 
the existence of normal extensions of partially known matrices based on earlier and some new results for 2x2, and 
2x3. We discuss also application of our criteriou to PPT entangled states of rank 5 in 2 x 4 systems. We relate these 
results to the theory of generalized concurrence [22| . 



Decompositions of density matrices. Physical states of composite quantum systems are represented by density 
matrices, i.e. Hermitian, positive definite linear operators of trace one, acting in the Hilbert space Tt = Ha ®Hb®..., 
which is a tensor product of Hilbert spaces corresponding to subsystems A, B, . . . of the considered system. In the 
following we shall be concerned with states of bipartite systems in a finite-dimensional Hilbert space i.e. described by 
positive definite Hermitian density matrices p = p^ > with Tip = 1, acting on the Hilbert space of the composite 
system Ti — Ha ®Hb- Without loosing generality we will assume that dim Ha = M > 2 and dim Hb — N > M, 
i.e. H — C (X) C N = C MxN . In the following we shall use the notation r{p) for the rank of the matrix p (which for 
a Hermitian matrix equals the number of its nonvanishing eigenvalues). 

Performing the spectral decomposition of p, 



where A; are (positive) eigenvalues of p, \ ipi) - its eigenvectors and r — r(p) - its rank. Defining | = \/A7| tpi) we 
thus decompose the nonnegative-definite Hermitian matrix p as a sum of rank-one operators, 



II. GRAM DECOMPOSITIONS OF DENSITY MATRICES 
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The decomposition of ((2|) into the sum a of rank-one operators is non-unique. Indeed the vectors 

r 

\® n ) = Y,W nl \V l ), n=l,...,K>r, (3) 
1=1 

lead to another one, 

K 

p = £l*»X*»l. (4) 

71=1 

involving K > r components, provided that the rectangular K x r matrix W fulfills W^W = I, where / is the r x r 
identity matrix, i.e. W belongs to the Vk,t manifold. In fact all possible decompositions Q of p into the sum of 
rank-one operators can be obtained from the spectral decomposition of p in such a way [H, [HI • 

Gram decompositions. Let {\E v )} v —i 
([2|) we obtain for the matrix elements of p: 

P„u = (E^\p\E„) - 



,MxN 



be a basis in 



iMxN 



Starting form the spectral decomposition 



i=i 



1=1 



(5) 



eC r , u£:=<*i|^>, (6) 



where (•, •) is the standard Hermitian scalar product in C K . It means that a Hermitian positive definite matrix is the 
Gram matrix (i.e. the matrix of scalar products) of the vectors w u defined above. A set of vectors {w„} fulfilling ([5]) 
we will call a Gram system for p, and we will say that it provides a Gram decomposition of p. If we do not insist 
that the vectors w v are elements of C we can construct other Gram decompositions of p, so in this sense the Gram 
decomposition is non-unique. Indeed, defining 

w' v = Vw v , (7) 

where w' £ C K , with K > r and V £ Mjf Xr - a rectangular matrix fulfilling V'V = /, we have p^ = (w'^, w' v ). It is 
easy to prove that all Gram decompositions of the matrix p are obtained by the transformation J7|) from the spectral 
one and ©. In particular two Gram systems calculated from two decompositions ^ and (0| are connected by 
the relation with V = W (cf. Q). In the following we will also use the fact that if two sets of vectors w' v and w", 
v = 1, . . . , r, wl S C K , w" £ C K are Gram systems for the same positive-definite matrix p, i.e. \w', u>l) = (w^, w ")> 
then there exists a unitary U acting in C K such that w" = Uw' v for v = 1, . . . , r. 

Gram decomposition in bipartite systems. Let us now take advantage of the fact that p acts on a tensor 
product space, i.e. we chose the basis {| £v)}„=i,....Af X jv in the form of product states | E v ) — \ e TO )®| f n ) =: | e m ®/„), 
rn = 1, . . . , M, n = 1, . . . , N and repeat the calculation of ([5]), 

r r 

Pij,mn = (e, ® fj\p\ e m ® fn) = ^2 (e, ® /j|*;)(*il e m ® /n) = X^^i'^™" = ( 8 ) 

/=i (=i 

where now w mn , m — 1, . . . , M, n = 1, . . . , N are vectors in C r , with components 

u>Ln = <*i|em<8/„}. (9) 

If we assume that p is of maximal rank r = MN (which for M > 3 we take for granted in the following), then w mn 
are linearly independent. In particular, for any fixed fh £ {1, . . . , M} (n £ {1, . . . , N}) the vectors »is n , n= 1, . . . , N 
(w m n, m — 1,...,M) form a set of N (M) linearly independent vectors in C K , respectively. In the special case 
M = 2, the first of the latter statements can be also assumed to hold, provided p is (non-trivially) supported in the 
2 x N space. If this statement was not true, the density matrix p would have a product vector in its kernel. In such 
situation either p is entangled and not PPT, or if it is PPT, than it can be represented as a mixture of a separable 
part and a density matrix supported in 2 x (jV — 1) dimensional space (for proofs and details see Ref. [2J|). 



w, 
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Relations between various Gram decompositions. Let {v n } n= i n form an arbitrary set of linearly inde- 
pendent vectors in C K (in particular, in the light of the above remark, we can choose v n = Wi n ). We can always find 
a family of linear maps F m : C K — > C K , m= 1, . . . , M, such that 

w mn = F m v n , (10) 

Note that such -F^'s are uniquely defined only on the iV-dimensional subspace of C K spanned by {w n }n=i,...Ar- If we 
apply the same procedure to the decomposition (j4]) we will arrive at 

w'mn = F L v 'm ( n ) 

with v' n , n = 1, . . . , N some linearly independent vectors in C K and appropriate F^ n . 
The vectors w mn and w' mn are connected by 0, i.e. 

F L V 'n = W 'mn = Vw mn = VF m V n . (12) 

Since both sets {v n } and {v' n }, n — 1, . . . , N are linearly independent in, respectively, C' r and C K , there exists a K xr 
matrix V of maximal rank, such that v' n = Vv n for n — 1, . . . , N . Consequently 

V^F' m Vv n = V*VF m v n = F m v n , (13) 
where we used V*V = I. We have thus shown the following 

Theorem 1 For the two decompositions of the Gram vectors of the form 110\) and 111]) stemming from two decompo- 
sitions of p into rank- one operators of the form 0) and there exist two K x r matrices V nad V , the former of 
rank r and the latter fulfilling V^V = I, such that on the space spanned by (arbitrary chosen) N linearly independent 
vectors v n € C r the equality 

V^F; n V = F m (14) 
holds. □ 
Obviously V depends on the choice of {v n }- 



III. SEPARABILITY PROBLEM 



Gram decompositions for separable states. Our goal is to characterize Gram decompositions for density 
matrices of bipartite separable quantum systems. Recall that separable defined on 'Ha®'Hb systems are characterized 
by the following 

Definition 1 A state p is separable if and only if 

k 

P = ^Pipf®pf- (15) 
i=i 

where ^2 { Pi — 1, Pi > 0, whereas pf and pf are states on TLa and Hb, respectively. 

The above expression means that p can be written as a convex combination of product states. 

To achieve the goal observe that performing a decomposition of the type ([2]) for all matrices pf and pf in (fT5]) 
and taking into account positivity of the coefficients pi we obtain that a state p is separable if and only if it can be 
decomposed in the form of K rank-one operators proportional to projections on simple tensors. 

K K 
k=l k=l 



where | <p k ) £ Ha, | ^fc) € H B - 
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Calculating matrix elements of p in local bases {| ei)}i=i,...,M and {| /<)}<=i jv in 7^ = C M and 7is 
respectively, we obtain 

Pij,mn — \ e i <8> /j |p| e m f n ) 
A 



!=1 
if 

E — /i 



ll / / t \ 



(17) 



where now 





' w' 1 ' 








w' = 


w' 2 


















w' K 









<p'l 





- w - 











(18) 



fo|e m >, i/>?:={tln\f n ). (19) 



Reformulation of the separability problem. From (fT8|) is is thus clear that for a separable state p on C M ®C N 
which can be decomposed into the sum of K rank-one product operators (fT6|) . there exist D m 6 M.kxK > to = 1, ... , M, 
D m - diagonal, and v„ € C x , n = 1, . . . , N, such that 

Pij,mn = (DiVpDmVn). (20) 

Eq. (j20|) is also a sufficient condition for separability, ie. if there exist v' n € C K , n = 1, . . . , N and diagonal D m 6 
M^xjf, to = 1, . . . , M such that ([20]) is fulfilled, then p can be decomposed into rank-one separable states (fT6|) with 

,1/ 



|W> = E (An)„|e m ), (21) 

m— 1 

TV 

1^) = EWI/«>- ( 22 ) 



i.e. p is separable. 
Indeed, define 



A' 



/0 = E|^i®^)<W®^l, (23) 



with if i and i/'z defined by (|2ip and (|22p . Then p is separable and an elementary calculation shows that pij, mn = pij, 
and thus p = p. 

Summarizing we can formulate thus the following theorem 

Theorem 2 A state p is separable if and only if there exists a Gram decomposition of p, 

Pij.mn — ij i ^mn ) ) 

Wij € C K for some K, such that 

ivy = DiVj, 

with N vectors {v±, . . . vn} € C K and M diagonal matrices D\, . . . Dm acting as operators on C K . 



Observe that we can assume that all diagonal matrices Di are nonsingular. Indeed from (fl8|) their diagonal elements 
are equal to the projections of the vectors \4>i), which constitute (a part of) the decomposition, onto the basis vectors 
| e m ). If any number of them vanish we can always adjust slightly the basis to make them taking non-zero values. 
Invoking now Theorem [T] we obtain a 
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Corrolary 1 A state p of the full rank r = MN is separable if and only if for some K > MN, there exist K x r 
matrices V and V of which V is of maximal rank and V'V = I, and diagonal K x K matrices D\, . . . Dm, such that 

V^D m V = F m 

holds on the space spanned by v n , n = 1, . . . , M , where w mn — F m v n is a Gram system {Pp for p calculated from its 
spectral decomposition 0). □ 

Before proceeding let us make a remark. Observe namely that in terms of the Gram decomposition (|20p of a 
separable p the operation of partial transposition in TLb 

p = J2p<p? ®pf^ p Tb = E^ ® ( 24 ) 

1=1 i=l 

ie. Pij,mn Pin,mj, corresponds to the complex conjugation of the frame {vi, . . . , i>n}, ie. 

v N } —>{«!,..., vn}, (25) 

whereas a similar operation performed in TIa 

k k 

p = Y,PiPt ® pf - ^ = E»(^) T ® ^ ( 26 ) 

i=l i=l 

16. Pij_rrm 1 * Pmj.im Consists m 

{Di, . . . , D M } — > {Si, • ■ ■ ,Sm}. (27) 

Indeed: 



^ n = P»n,mj = (DiV n ,D m Vj) = (v n ,DiD m Vj) = (DiD m Vj,V n ) = (DiD m Vj,V^} 

= (D m DiVj,v^) = {DiVj,D m v^), 
pj 3 \mn = Pmj.in = {D m v 0l DiV n ) = (vj,D m DiV n ) = (vj,DiD m v„) = (DiVj,D m v n ), 

where we used the fact that diagonal matrices commute and their Hermitian conjugation reduces to the complex one. 

In Appendix [X] we discuss in detail the results of this section applied to an example of two qubit states, the so 
called Werner states. 



IV. SEPARABILITY AND FULL FAMILIES OF COMMUTING NORMAL MATRICES 

K -separability and FFCNM. In forthcoming publications we will present applications of Corollary [1] to char- 
acterization of the bipartite entanglement for arbitrary systems. In the present paper we would like to concentrate on 
the separability of systems with M = 2, but before that we would like to use our results from the previous sections to 
formulate a novel SNC for separability. Let us assume that the investigated state p is ^-separable, i.e. there exists a 
decomposition into exactly K rank-one product operators (flU]) and p can be cast in the from (f2"0")) . Since we assumed 
that p is of maximal rank r = NM we have necessarily K > MN. From the previous consideration we now that the 
Gram vectors calculated with the help of this decomposition have the form w' mn — D m v^. 

For an arbitrary decomposition of p into exactly K states given by (0J (where | $/) need not to be product states), 
we obtain another Gram system w^ n = ($z|e m ® /„) € C K . 

The vectors w' mn and w!^ n as forming two Gram systems for the same matrix p are connected via a unitary 
transformation 

«d = u w Ln = UD m v' n . (28) 
Taking the above equality for two pairs of indices (to, n) and [k, n) we obtain: 

M mk w£ n = w£ n , (29) 

where 

M mk ^VD m {D k Y x \J^. (30) 

Remember that without loosing generality we can assume nonsingularity of all matrices D n . Consequently M km are 
also nonsingular. 
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M km ,Mt 



km 



0, and mutually com- 



5 ATC for separability and FFCNM. 

muting, [Mf. m , Mi m i] = 0. Both observations can be easily proved using the facts that all matrices D rn are diagonal 
and U is unitary. The above reasoning is summarized in the form of the following 

Theorem 3 A necessary and sufficient condition for K -separability of p is the existence, for an arbitrary decompo- 
sition Q), of a full family of M{M — l)/2 normal, commuting K x K matrices M km fulfilling i29\) where w' mn are 
appropriate Gram vectors for the decomposition \29\). 



Necessity of the condition follows from the above remarks, and to prove the sufficiency let us assume that (129)) is 
fulfilled for some family of normal, commuting matrices M km . It is a standard fact from the linear algebra [23] that 
all matrices in such a family can be simultaneously diagonalized by a single unitary transformation, 

uHd km U = D km . (31) 

According to the previous remarks we assume that M km and, consequently, D km are nonsingular. Now from (|29[) and 

(ED 



D km U^w^ n = UHl km UU^w^ n = U^M km w^ n = rfw^, (32) 

and defining v n := Ww" n , w kn :— U^u>l! n , we obtain 

w kn = (-Difc)" 1 "™- (33) 

The vectors w kn are Gram vectors for p as they are obtained by a single unitary transformation form the vectors 
w' k ' n constituting some Gram decomposition of p. Equation (|32l) reveals their structure in the form sufficient for the 
separability of p according to Theorem [2] 

V. SEPARABILITY IN 2 x N SYSTEMS AND NORMAL EXTENSIONS 

The Theorem 3 simplifies significantly for 2 x N systems, because FFCNM consists of a single matrix, which has 
to fulfill 

M< = <. (34) 

In the following we will use 0, 1 instead of 1, 2 for numbering the components on the qubit side, which is more in accord 
with the custom to denote the basis states by | 0) and 1 1). From here we do not need to assume the nonsingularity 
of p - see the remarks preceding the formula (flOl) . 

In this section we study the consequences of ()34|) . On one hand we use the present formulation to obtain particularly 
simple proofs of known separability criteria. On the other hand, we use known separability criteria to obtain non- 
trivial statements concerning existence of normal extensions of matrices, whose matrix elements are only partially 
known. 

Canonical forms and PPT condition. Let us consider p in the canonical form [24| (see also Appendix IB"]) 



A B 
fit i 



(35) 



where the positivity of p implies A — BB^ + AA^, where A is some N x p matrix, with p > r(AAT). Obviously, 
p = 1 necessarily, when r(AAT) = 1; also one can always take the minimal p = r(AAT) = \ We represent AA^ = 
ELi|A„)(A n |. 

Similar considerations concern the partially transposed matrix, which reads 



p Ta 



A St 
B I 



(36) 



We consider here only the nontrivial case of states with the positive partial transpose (PPT states) - states which are 
not PPT are not separable. The positivity of p TA requires now that A = B^B + A'A, where A is now a p x N matrix, 
withp > r(A' l 'A), and having analogous properties asp introduced above. We represent AT A = J2 n =iP\ ^n)(A n |. The 
PPT condition can be thus stated as A — B^B > 0. More precisely, it must hold 

A = BB f + AA f = B^B + A f A, (37) 

which implies that given A, A are not independent, and related by the above constraint. 
Let us now discuss several examples to show how the novel entanglement SNC works. 
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Rank N matrices. The results of Ref. [2J| indicate that rank N PPT states are TV-separable. The matrix M = B 
then, and [B, B f ] =0. 

The case p = p TA . From Ref. [24| we gather also that when p — p TA , then p is 27V-separable. In this case 
B = B\ A = At and the matrix M can be written as 



M = 



B A 
At s 



(38) 



with N x N matrix S to be determined. Obviously, taking s Hermitian provides the desired normal extension of M. 

The case of 2 x 2 and 2x3 systems. In the two qubit, or qubit-qutrit case, every separable matrix is in- 
separable, where K = m.ax(r(p),r(p Tjl )). In particular for the full rank r(p) = 4 (r(p) = 6), K = 4 [2(| or K = 6 (as 
shown in Appendix B), respectively. We have then 

Corrolary 2 For N — 2,3, an arbitrary N x N matrix B, and an arbitrary p x N matrix A constrained by fr37\) , the 
matrix 



M 



B A 
A s 



(39) 



a has a normal extension, i.e. there exist a p x p matrix s, and a N x p matrix A constrained by \31\ )such that M is 
normal. This holds in particular for minimal p = min 0, r(p) — N. 

Edge PPT entangled states for N = 4. Per hap s the most interesting are applications for PPT entangled 
states, and in particular for the so called edge states [24|, i.e. PPT states that cannot be represented as a mixture of 
a PPT state and a separable states (no separable part can be subtracted) . Such states are extreme examples of states 
to which the range criterion of P. Horodecki [27j applies. For N = 4 such states may have rank 5, or 6 (and similarly 
their partial transpose). From the analysis of the Appendix B we infer that if p is a separable state of rank 5 such 
that its partial transpose has rank 5 (6), then it is 5-separable (6-separable). In the case r(p) = r(p TA ) = 5, both A 
and A have rank 1; we denote A by | A), and At by | A). We get then 

Corrolary 3 A PPT state p such that it and its partial transpose have rank 5 is separable, iff there exist a complex 
number s such that the matrix 



M 



B | A) 
<A| 8 



(40) 



is normal, which assuming that \37\ ) holds, requires that 

(B- s )|A) = (5t 



( )|A). 



This condition is equivalent to the range criterion. For the particular example p<jj of the 2x4 state analyzed in the 
seminal 1997 paper [23], it is very easy to analyse, as we show in Appendix C. 

This analysis may be extended to the rank 5 states, with the partial transpose of rank 6, for which Eq. (f3"T)) becomes 



BB* + | A) (At | = B*B + | A\) (Ai | + | A 2 ) (A 2 |. 



(41) 



We have in this case 



Corrolary 4 A PPT state p of rank 5, such that its partial transpose has rank 6, is separable, iff there exist complex 
numbers a, (3, such that \a\ 2 + |/3| 2 = 1, and a 2 x 2 matrix s such that the matrix 



M 



B a\A) [3\ A) 

(Ai| Sn Si2 
(A 2 | S 2 1 S22 



(42) 



is normal, which assuming that |^j[ ) holds, requires that 

(B-sii)|Ai)-«2i|A a ) 



(a£t 



(B - s 22 )\ A 2 ) - s 12 | Ai) = (fiB* - as* 21 - (3s* 22 )\ A). 
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Before we end this section, we would like to stress that obviously the above discussion of normal extension of M 
applies also to M x N systems, if we focus on a single relation of the type ([29"]) , such as say 

M w w^ n = <. (43) 

The analysis pertains then to the study of separability on a particular 2 x N subspace of the full Hilbert space. In 
this sense it is somewhat similar to the theory of generalized concurrences of Ref . [22| • 



VI. SUMMARY 



We have presented a new approach to the separability problem by reformulating it in terms of existence of separable 
Gram decompositions of density matrices in auxiliary space. The existence of such Gram decompositions is equivalent 
to the existence of a full family of commuting normal matrices that relate components of Gram vectors. We have 
presented many examples and applications of this method mainly to the 2 x N systems. Several known separability 
criteria can be, on one hand, reproduced with this method in a particulary simple way, and on the other, can be used 
to derive nontrivial statements about the existence of FFCNM. 
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APPENDIX A: WERNER MATRICES FOR TWO QUBITS 



As an example illustrating the results of Section HTT1 let us consider a one-parameter family of states for N = M = 2 

1+p 2p ■ 
j 1-p 
1-p 
2p 1+p. 



(Al) 



the so called Werner states [J]. The parameter p takes the values from the interval [0, 1]. One finds easily the spectral 
decomposition p = XT=i l^/X^il w ^ tn 



*!> = 



and calculates the Gram vectors © 



























V 8 





fi^p 







V 8 


,l* 2 ) = 


. 1 *3> = 


V 8 


. 1 *4> = 


/i-p 





n^p 





V 8 





V 8 




V 8 













1+3 P 



W12 = 



s 



s 





,W 2 1 = 



We chose v\ — wn and v-i — W22 which allows to take 



Ft 



10 
10 
10 
1 






-1 


-1 
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, W22 = 



1+3 p 



l+3p 
1-P 



1-P 
l+3p 





(A2) 



(A3) 



(A4) 
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Only when p < 1/3 the state p is separable. For these values of p one finds an explicit Gram decomposition (|20| of 
p with 



10 
10 
10 
1 



and 



£>2 



(i-i)(yr^+VFFT) 
%/2(yr^+y3^+T) 







v / 2(v / T Tr P-v / 3p+T) 



(yi^+V5+i) ™ 

t= -e 2 

4V2 

(v/T^-V^P+T) ai 

-e 2 

(v^+VSFFT) _ni 

?= -e 2 

4y2 



-e 2 



(i-i)(yT^3F+VF+T) 

\/2(yT=^+y3p+T)~ 

(!+»)( yp+T-yT^) 
V2(yT=£-V3p+T) ' 



(yp+T-yT=3F) 2s* 
7= -e * 

(yT=3p+yp+T) _3 ZL i 

- —p 4 

4y2 

(Vp+T-\/T^3F) 3^ 

- t= -e 4 

4\ /2 

(yT=3p+yp+T) 

-p 4 

4v^ 



For the particular choice of F\ , F2 , v\ , and V2 the matrices V and 1/ (cf . Corollary [T]) are given as 



V = V 



/T+p+Wl-3p 



yi — 3p — iyi+p 



-yi-3p+i:yT+j; 
ylyT^p 



/T+p+iVl-3p 



yi— 3p-^yi+p 
ysyT^p 



-yr+P+jyi-3p 
ysyr^ 

for which one easily checks D m V = F m on span(vi, V2) 



1 — y i—3 p+iyi+p 

2 ysyr^p 



(A5) 



APPENDIX B: A SHORT GUIDE TO 2 x N SYSTEMS 



From Theorem [3] it is clear that investigations of separability can be simplified if we know a priori the order 
of separability K of the given state (i.e. we know that it is if-separable) . We do not have any general tool for 
determining exactly the order of separability for arbitrary separable states before finding their actual decomposition 
into pure products (and even if we find one, to establish the order of separability we still have to prove that the found 
decomposition involves the minimal number of components). Here we present some exact results concerning orders 
of separability in the case of 2 x TV systems for low values of N. 

Canonical forms. Let p be an arbitrary density matrix of a bipartite 2 x N system 



P 



A B 

St C 



(Bl) 



where A, B and C are N x N matrices, A and C are hermitian due to hermiticity of p. Positive definiteness of p 
implies A > 0, C > 0, and A - BC~ X B^ > 0. The matrix C is nonsingular since p, by assumption, is of maximal 
rank. By an invertible transformation 



(i®c- 1 ' 2 ) p (i®c- 1 ' 2 ') 



(B2) 
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we bring p to the canonical form [24 



P 



A B 
St I 



(B3) 



For simplicity of notation we kept the same symbols A and B to denote the appropriate blocks of p despite the fact 
that the original blocks defined in (|Bip are altered by the transformation (|B2|) . Such a transformation changes, in 
principle, the trace of p, but since the normalization of the trace does not influence separability properties we will 
wave this point aside. The positivity conditions of p reduce now to 



A = BB^ + AA t , 



(B4) 



where A is some p x N matrix. 

A necessary criterion of separability is the non-negative definitcncss of the partially transposed matrix, which for 
the case of a 2 x N system (2N) is defined as 



P 



A £t 
B I 



(B5) 



From now on we will assume thus that both p and p TA are positive-definite (otherwise p is not separable). Positive 
definiteness of the partial transpose of (|B3[) demands 



A = B^B + A^A, 



(B6) 



for some A. 

Decompositions of p. For our purposes it important to consider particular decompositions of p for K = N + p, 
and construct the "known" part of the matrix M that fulfills Eq. (j2"5|) . From the canonical form and Eq. (|B4|> . with 
A = (| Ax),..., | Ap)), it is easy to see that p = J2k=i l$fc)( $ fc|, with 



**) = \0)®\k) + \l)®B\k), 



for k = 1, ... , N, and 



1)®I A 



k-N), 



for k = N + 1, . . . , N + p. From this particular form we read the components of the vectors Wq u , w" n 



for k — 1, . . . , N, and zero otherwise. Similarly 



for k = 1. 



, N, and 



«J fe = B nk , 



«J fc = (n\A„- N ), 



(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 



for k = N + 1, . . . , N + p. Obviously, for the particularly simple form of uiQ n , Eq. (|29|) . determines only the first N 
columns of the matrix M T , which are 



M T = 



B T ? 
At ? 



(B12) 



where the other entries are at this moment not known. Transposing and using the PPT constraint (I37|) . we indeed 
obtain that 



M 



B A 
A S 



(B13) 



where for given B and A, the matrix A is constrained only by the condition (|57| . and S is completely arbitrary. This 
form of M is intensively used by us in the section on separability in 2 x N systems. 
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Edge states. For the purpose of this paper we remind the reader the basic concept associated with the edge 
states. First, we remind pit \2&L I2H] that if | e, /) (or any vector, in fact) is in the range of p, then we can write 

p = p' + \\ej)(ej\, 

where p' > provided A < l/(e, /|p _1 |e, /). When the equality holds, the rank of p' is smaller that the rank of p by 
1. 

From this observation follows 

Corrolary 5 A PPT state p is an edge state if there exist no product vector | e, /) in its range, such that \ e* , f) is 
in the range of p TA . 

Subtracting projectors on product vectors as in Refs. 0, [H| allows to determine the minimal number of projectors 
on product states necessary to decompose a separable state. We remind the reader: 

The case 2x2. In this case K = max(r(p), r(p TA )). This result stems from (26[. In the following we shall use 
notation (p, q) for the case of r(p) = p, r(p TA ) = q. Let us consider the case of full ranks (4,4). First we show that 
we can find the product vector for which A = l/(e, /|p _1 |e, /) = l/(e*, /Kp^)" 1 ^*, /), so that subtracting projector 
on this vector reduces ranks to (3,3). To this aim we suppose p — J2k=i l e fc> fk)( e k, fk\ is separable, and that for all 
|e,/), it holds (e, f\p~ 1 \e, /) < (e*, f\{p TA )~ 1 \e* . f). Inserting into this inequality \ek,fk) and summing over k, we 
get a contradiction Tr{I) = 4 < 4 = Tr(I). In the same manner we prove that the opposite inequality can not be 
fulfilled by all product vectors. Thus either all product vectors fulfill the equality, or there are at least two product 
vectors for which the inequality takes opposite signs. But then from the Darboux property and the fact the the set of 
all product states is connected, we gather that there exist a product vector for which the equality holds. In the next 
step we reduce one rank to 2; this, however implies that so does the other rank, since from general theory of Ref. [24| 
it follows that rank N PPT matrix in 2 x N systems is ./V-separable. 

The case 2x3. This problem was partially addressed in the thesis of G. Vidal 29]. The proof here is new. We 
start with the full ranks and using the same argument as above we reduce the ranks to (5,5). The argument may be 
then repeated but with a certain care. Now we suppose that p = Y]u—i \ek, fk)(ek, fk\ is separable, and that for all 
I e, /) in its range, and such that | e*, /) is in the range of p TA , it holds (e, f\p~ 1 \e, f) < (e*, ,/|(p TA ) -1 |e*, /). Again, 
inserting into this inequality | €k,fk) an d summing over k, we get a contradiction Tr{I R ^) = 5 < 5 = Tr{I R ^ p T A ^), 
where Ir(p) denotes identity on the range. We may again evoke the Darboux property, but to this aim we need to 
prove that the set of product vectors on question os connected. Let $ be a vector from the kernel of p and $ from 
the kernel of p TA . The product vectors we look for have to fulfill (\& |e, /) = 0, (3>|e, /) = 0. These equations can be 
regarded as two linear equations for a three-component vector | /), parametrized by the vector | e) — | 0) +a| 1), which 
we have parametrized by complex number a in some basis. Obviosly, | /) is a unique function of a and by scanning a 
over the complex plane we can reach any of these vectors in a continuous way. Darboux theorem says then that there 
exist a product vector for which equality holds (e, f\p~ 1 \e, f) = (e* , f\(p TA )~ 1 \e* , /), and we can reduce the ranks to 
(4,4). The next step is as above: reduction of one of the ranks to 3, implies the same reduction for the other. The 
reason for that is that all rank 3 states in 2 x 3 systems are 3-separable. 

The case 2x4. It is also possible to determine what is the minimal number of terms in the separable decomposition 
for the states of low ranks. In this paper we consider two cases: (5,5) and (5,6). In the (5,5) case there are three 
vectors | in the kernel of p, and another three vectors | <&i) in the kernel of p TA . We look for | ek, fk) such that 
fi>i\e,f) = 0, ($j|e*,/) = for all i — 1,2,3. These can be regarded as six linear equations for a four-component 
vector I /) . They have solutions provided three 4x4 determinants (constructed from the first three and one of the 
last three equations) vanish. These determinant constitute three polynomials of 3rd order in a and first order in 
a*. Eliminating a* from them we obtain that two polynomials of 6th order in a must vanish. Subtracting them 
with appropriate coefficients, we conclude that a polynomial of 5th order in a must vanish, i.e. there are at most five 
product vectors having the desired properties. This implies that if p is separable, then it is 5-separable. 

Similar analysis can be done for the case of the state p with the ranks (5,6). We end up then with one polynomial 
of 6th order in a, i.e. we have at most six solutions, ergo if p is separable, then it is 6-separable. Note, that the states 
with ranks (5,6) are either separable, or entangled edge states, or mixtures of ranks (5,5) edge state with a single 
projector on a product vector from the range of p. 

Unfortunately, only upper bounds on the number of product states in an decomposition of separable states are 
known for p's of higher ranks. In particular, Caratheodory theorem (for proof see [271 ]) gives a general bound equal to 
the square of the dimension of the Hilbert space, i.e. in the present case (2 x 4) 2 = 64, implying the every separable 
state is 64-separable. 

For the states with ranks (5,7) it can be shown that there exists in the range of p a product vector (| 0) + a\ 1))| /), 
such that (| 0) + a*\ 1))| /) is in the range of p TA . It is easy to see that the condition that these product vectors 
are orthogonal to the corresponding kernels of p and p TA , leads to 4 linear equations for 4 components of | /) . The 
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solutions of such equations exits if the determinant of the corresponding matrix vanishes. This matrix has three rows 
linear in a and one row linear in a*, so that the determinant equations has the form 

W 3 (a) + a*V 3 {a) = 0, (B14) 

where Ws(.) and V 3 (.) are polynomials of third order. Let us replace a — > rs, a* — > r/s with r > and s complex, 
and treat Eq. (|B14[) as an equation for s(r) (i.e. treating s as parametrically dependent on r), 

sW 3 (rs) + rV 3 (rs) = 0. 

We will show that this equation has at least one root a = rs with |s| = 1, i.e. with a* — r/s. To this aim we consider 
the asymptotic behavior at r — ► oo. It is easy to show that the above equation has three roots Si = 0(l/r) — ► 0, 
i = 1,2,3 and one root S4 = O(r) — * 00. Analogously, for r — > it is easy to show that the equation has three roots 
Sj = 0(1 /r) — * 00, i = 1, 2, 3 and one root S4 = 0(r) — > 0. All that implies that when we continuously change r from 
to 00, one the the three "large" roots must become "small". From continuity (i.e. again from the Darboux property) 
we get that for some r = rg, the |s(ro)| = 1. Unfortunately, we cannot say much more about the total number of such 
roots. Solving Eq. (|B14[) with respect to a*, complex conjugating the result, and stacking it back into Eq. (|B14[) . we 
obtain an equation for a* of 10th order, which indicates that there are not more of than ten roots of Eq. (|B14I) . 



APPENDIX C: HORODECKIS 2x4 EDGE STATE 



In this Appendix we show how our method work for the famous state P97 introduced by P. Horodecki in the seminal 
paper [27]. In our notation this state has 



B = 






1 





0" 








1 











1 


















(CI) 



and 



A) 



V(l-6)/26 



y/(l + b)/2b 



A) 



y/(l + b)/2b 



y/0--b)/2b 



(C2) 



This state has rank 5 (as its partial transpose) and is an example of an edge state |24j . As pointed in Ref . [27J there 
exist a unitary matrix K, such that K 2 = I, KBK = B>, and K\ A) = | A). The condition of existence of the normal 
extension from Section VI reads then [B — s)\ A) = K(B — s*)\ A), i.e. 



-sy/(l + b)/2b 


V( l - b)/2b 
-<V(l-6)/26 



-s* y (1 - b)/2 b 




(C3) 



which has only the two solutions s = 0, b — 1, and the limiting case b 
exactly the two instances in which the Horodecki state is separable. 



s*^/{l + b)/2b_ 

0, with an arbitrary real s = s*. These are 
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